Role of Bound Magnon in Magnetic Domain Wall Motion 
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We report on a quantum description of the domain wall (DW) motion under a spin current. A 
bound magnon, which is the zero mode of D W, is found to play a dominant role in D W dynamics. 
The bound magnon acquires its inertia by the hard axis anisotropy and is a free particle even under 
the spin current. The full transfer of spin angular momentum from the spin current to DW via the 
bound magnon leads to the D W motion with the adiabatic velocity, decoupling of spin waves from 
DW, and no Doppler shift in spin waves. 

PACS numbers: 75.45.+j, 75.60.Ch 



Magnetic DW motion [l], 0] has attracted much inter- 
est to both experiments and theories, owing to its po- 
tential device applications such as DW logic [3j and DW 
memory [4J. DW motion under magnetic fields is now 
well documented and understood [l| in terms of classi- 
cal Landau-Lifshitz- Gilbert (LLG) equations. Recently 
many experimental groups 0, H, 0, @, S E3| have ob- 
served the magnetic DW motion under spin currents, 
theoretical description 
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However, theoretical description [11|, 1121 H 
very controversial over several issues. Even for a perfect 
ferromagnetic (FM) nanowire, the hard axis anisotropy 
was claimed [12| to induce the intrinsic pinning so that 
a DW does not move until a finite spin current is ap- 
plied. But this intrinsic pinnin g is not consistent with 
translational symmetry of DW [l4j. In the generalized 
LLG equation approach to the DW motion, the spin cur- 
rent gives rise to the so-called nonadiabatic 13. 14, lB| 
as well as adiabatic torques. The relative magnitude of 
nonadiabatic and Gilbert damping torques remains an 
unresolved problem. 

In this paper we study the DW motion under spin 
currents, using the one dimensional s-d model Hamilto- 
nian and the full quantum mechanical description of the 
DW motion. In the LLG equation approach, macroscopic 
magnetization is treated as a classical vector with fixed 
magnitude and two Euler angles. In this work we treat 
both electron and local spin systems quantum mechan- 
ically in order to address a microscopic mechanism for 
spin transfer from a spin current to a DW. We find that 
a bound magnon (with zero energy) in the DW plays a 
dominant role in absorbing angular momentum from the 
spin current and thereby in the DW motion. This bound 
magnon is the zero mode of magnetic soliton (DW) [lij ]. 
The mass of the bound magnon is derived quantum me- 
chanically for the first time and is in agreement with the 
classical DW or Doring mass [17J . Our study shows that 
the bound magnon is a free particle even under a spin 
current and a DW can start to move under any finite 
spin currents for a perfect FM nanowire. The bound 
magnon-mediated spin transfer mechanism is compatible 



with the balance of the nonadiabatic and Gilbert damp- 
ing torques. The role of the bound magnon is very similar 
to that of zero phonon mode in solitons of charge density 
wave state fljj . 

We consider the one-dimensional s-d model to study 
the DW motion which is driven by spin currents. Our 
model system consists of three parts: H = H e + Hs+H e s 
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(1) 
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H S = -jJ2Si-S i+1 -Aj2(Si-zf + Kj2(Si-yp) 

i i i 

H e S = —Jll'yi^ci " (3) 



Conduction s electrons are described by H e and are spin 
polarized by the Hund coupling H e s to the ordered local 
spins. Hs describes a system of local spins with an easy 
z-axis {A > 0) along the wire direction and a hard y- 
axis (K > 0). Ferromagnetically (J > 0) coupled local 
spins are assumed to have a transverse DW. DWs can 
be induced in the FM wire by ingenious experimental 
techniques. The mutually orthogonal unit vectors, x, y 
and z, define the laboratory frame. 

Representing the local spins in terms of two Euler angle 
fields, 9 and </>, the DW structure can be derived by the 
energy minimization. The ordered local spins lie on the 
easy x-z plane (<f> w = 0, 7r) and are rotated away from the 



easy axis by the angle 8 
DW position, A 



(z) = 2 cot' 1 e 



-(z-q)/A 



q is the 



y/Ja 2 /2A is the DW width and a is 
the lattice spacing between two neighboring spins. The 
rotation angle W at each spin defines the quantization 
axis along which local spins are aligned. 

Local spins may well fluctuate away from the ordered 
DW state. Small fluctuating spin fields can be repre- 
sented by small fluctuating Euler angle fields. But in- 
stead, we adopt in this work the small fluctuating trans- 
verse spin fields, Si X and Si y , which are a more natural 
description of spin fluctuations. Using the local coordi- 
nate frames defined by the local spin quantization axis, 
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the spin fields can be represented in terms of local trans- 
verse spin fields, 

x ■ Si = Si X cos w (zi) + S iz sin 6 w (zi), (4) 
y ' Si — Siy, (5) 
Z'Si = SizCosd w (zi) — Si X sm8 w (zi). (6) 



wave operators) and b w = J2i V'bC 2 *)^ ( a bound magnon 
operator). Since ipB and ipk's form a complete set of or- 
thonormal wave functions for Hd W , b w and b^s exhaust 
all possible normal modes of H and the inverse relation 
can be readily written down, 



h = 



Here S iz = ^JS(S + 1) - Sf x - Sf y . Expanding the 

transverse spin fluctuations away from the quantization 
axis, we find the magnon Hamiltonian H rnag = Hq + Hi, 
keeping only up to quadratic terms in transverse spin 
components. 



k 



(12) 



"o = EE 



i a—x,y 



2 {Sia Sl+la) 



A(cos 2 9 Z - sin 2 Oi)Sl 



In terms of normal modes b w and 6fc's, Hq is already 
diagonalized: H = J2k e k"kbk- Due to its zero energy, 
the bound magnon does not show up formally in Hq. 

Including the hard axis anisotropy (HAA) , the magnon 
Hamiltonian can be diagonalized as H mag — Hq + H sw . 
The bound magnon still remains as a normal mode with 
nergy E B = 0. 



and Hi — KJ2i where 9i — 9 w (zi). The transverse 
spins in the continuum limit (Si X — Si+i x ~ —adS x /dz) 
satisfy the following equations of motion, 



hd t S x = HdwSy + 2KS Sy, 

hdtSy — —HdwS x - 



(8) 
(9) 



Here H dw — — JSa 2 dl + 2AS(cos 2 9 W — sin 9 W ) and acts 
as a Hamiltonian for normal modes of magnons. 

Winter [ij} found two types of normal modes for the 
transverse DW: the magnon bound to the DW and the 
extended spin waves. The static Schrodinger equation of 
Hdw accommodates one bound state with energy e = 
and extended states of wave number k with eigen energy 
e k = JSa 2 k 2 +2 AS [ljj. 



ikA — tanh 



Vl + fc 2 A 2 



A e ik(z-q) 



(10) 

(11) 



iPb(z) is the normal mode with the energy eigenvalue 
e = 0. This normal mode is bound within the potential 
well formed in the DW. Due to localization in the DW, 
this magnon mode may be called as the bound magnon. 
The wave function ipg of the bound magnon is related 
to the DW structure function 9 W by its spatial derivative 
d z 9 w , which implies a translation of the DW. This bound 
magnon is none other than the zero mode [l6| of a DW, 
and tends to restore the translation symmetry of the FM 
nanowire. On the other hand, extended spin waves ipk(z) 
of wave number k (L: length of FM wire) are plane waves 
with reduced amplitude in the DW, are characterized by 
the excitation energy gap and deforms the DW. 

Local spin fields can be represented with the Holstein- 



Primakoff magnons: Si 
bt 



2S-b% bt 



2Sh, 



2Sb] and S iz = S - b\b t . With 



Si- = b\yj2S-b\bi'. 
identification of the normal modes, we can define the cor- 
responding magnon operators: bk = J2i^k( z i)bi (spin 



H B 

S 3 
S L 



■\KS(b w -blf . 



(13) 
(14) 



Note that the bound magnon has no spin component 
along the hard axis, S y = 0, but instead, its spin lies 
on the easy plane. S x has the zero mode, which means 
that spins can rotate freely on the easy plane so that the 
DW can be shifted freely along the FM wire direction. 
No zero mode in S y simply reflects no free rotation of 
spins away from the easy plane. On the other hand, 
the spin waves have spin excitations along two trans- 
verse directions with the increased excitation energy gap 
E k = yje k (t k + 2KS) under the HAA. 



H sw = E k a k a k , 



oc ipk{z) 



u x (k) 

Uy{k) 



(15) 



(16) 



Here is the spin wave boson operator under the HAA 
and a linear combination of b k and 6^ fe , and u x / y corre- 
sponds to the amplitudes of transverse spins. 

For the ordered local spins, the spin texture can be 
described by the magnetization unit vectors, mi's, where 
m, = fn(zi) and fh = z cos 9 + x sin 9 cos (f> + y sin 9 sin tfi. 

Under unitary transformation Ui — exp (— ^9i^>i • a J , 
which rotates the quantization axis of conduction elec- 
trons at site i from the z axis into fhi (ci a i— » di a = 
u*p a Ci@), i.e., u\a ■ rhiUi = a z , the Hund coupling H e s 
is diagonalized and results in spin polarized conduction 
bands. The kinetic term in H e introduces the current- 
spin coupling H c s or the Berry phase term fill]. 



H c s = frv s S E(! ~ cos 0*)Qg0t- 



(17) 



Here v s = and I s is the spin current flowing in the 

system under electric field. I s is computed from ther- 
mal average of I si = || ^ a a[d\ +la d la - d\ a d l+ i a ], which 
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measures the spin polarized electric current from i to 
H c s can be written in a compact form as 



H cS = v s V, 



(18) 
(19) 



V is the DW linear momentum [21( or the generator of 
DW translation as will be shown below. Here angles are 
field variables. 

We now prove that V is the generator of DW trans- 
lation or the linear momentum for DW. For this pur- 
pose we consider the DW spin texture |^({zi}) >= 



n \Si 



>= Yli U i{mi)\ s uZ >! where U t i 



cxp (—i9j(f>j ■ Sij rotates the orientation of Si from the 

z axis into fhj. The DW state shifted to right by a 
lattice constant a can be written as ^({zi — a}) >= 
U.i\Si',Ai-i >= UiUiim-iWlim^iizi}) >. Writ- 
ing |*({z 4 -a}) >= exp(-ifV) |*({^}) >, we can iden- 
tify V as 

V = h Si ■ [rhi(l — cos6i)d z 4>i — fhi x d z rhi] (20) 

i 

This quantum definition of V can also be obtained from 
Eq. (fl9|) by allowing small fluctuating angle or spin fields 
as in the normal mode expansion. Angles or m; represent 
the DW solution. The first term is c-number [S% ■ fhi — 
Si Z ~ S), while the second is the quantum correction P 
and 



P = -hJ2d z 9 t S iy , 



(21) 



for the transverse DW in our case. 

Under the finite spin current, w {z) now becomes dy- 
namical and thus, the DW position q is time-dependent. 
Two coupled Eqs. © and ©, under the current-spin 
coupling (fPB")) . are modified by two effects: dynamic 6 W 
and the spin current. In the rotating frame about the 
hard axis or the y axis, the equation of motion for an op- 
erator A is ihdtA — [A, H c g], where the effective Hamil- 
tonian H c ff = H + Hi + H c $ 
contribution from rotating angle 9 



h6Sy has an additional 



hd t s x = H dw s y + 2KS s y -hS(v s d z e + d t e), (22) 



hdfSy — —Hd w S x 



(23) 



In general, the spin waves are coupled to the DW. Nor- 
mal modes under spin current are decoupled from DW 
only when t| = v s , i.e., the DW moves with the adia- 
batic velocity v s . If the DW absorbs with full efficiency 
the spin angular momentum transferred from the s pin 
current, there will be no Doppler shift 11, HH, [3, 0] 
in the spin wave energy spectrum. If not, spin angular 
momentum from the spin current will be transferred to 
exciting spin waves. 



The DW dynamics is determined by the bound magnon 
Hamiltonian, Hg and H c s- Since d z cx tJjb is finite only 
near the DW, the main contribution to P comes from 
spins in the DW and the number of contributing spins 
is roughly A /a. Furthermore, owing to Hdw sin 9 W — 
0, we have the identity [P, H mag \ — such that P is 
the constant of motion. This is a simple mathematical 
manifestation of translational symmetry for a DW in an 
infinite FM nanowire. P can be represented in terms of 
the bound magnon as 



P = -ih\ 



b w ). 



(24) 



Hb can be interpreted as the kinetic Hamiltonian of the 
bound magnon by noting that Hb can be written in terms 
of P as 



Hr 



KAa 



-P 2 = 



P 2 



2M dv 



(25) 



A" A a 



where the bound magnon mass is defined as Mdv 
The DW or Doring mass defined in the classical approach 
171 ] is none other than the mass of the bound magnon or 
the zero mode in a ferromagnetic DW. 

The effect of the hard axis anisotropy K is threefold. 
For extended spin waves, their energy gap is enhanced 
such that they become much harder to excite. The hard 
axis anisotropy confines the bound magnon to have spin 
components only on the easy plane, but no component 
along the hard axis. The bound magnon acquires its 
inertia due to the hard axis anisotropy. 

Dropping c-number from V, the current-spin coupling 
Eq. (fl8|) becomes 



H 



cS 



= v,P, 



(26) 



which is the same for both <f> w — 0, 7T. Note that P(4> w = 
7r) = —P{4>w = 0). Under spin currents, the system 
retains a translational symmetry. The action of H c $ on 
the DW can be most easily understood in terms of the 
Schrodinger equation, 

ihd t \^({zi},t)> = (H B +H c s)\V({zi},t)> .(27) 

Denoting the DW state as |\I/({zi}) > when v s = 0, 
we find that >= exp(—iv„tP/h) \V({zi}) >= 

\^({zi — v s t}) >. The DW motion with velocity v s is 
induced by the spin current. 

The intrinsic pinning was claimed [12j to be induced 
by HAA of a perfect FM nanowire. Spin currents ro- 
tate the local spins away from the easy plane and the 
HAA field acts as a blockade 12] to the DW motion and 
generates the intrinsic pinning. The quantum approach 
clearly shows the absence of the intrinsic pinning and is 
consistent with the translation symmetry of DWs under 
spin currents. The DW absorbs the spin angular momen- 
tum from the spin current via the bound magnon, and 
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thereby avoids the tilting of local spins away from the 
easy plane. 

There is (no) translation symmetry for DW in FM 
wires under spin currents (magnetic fields). The energy 
damping torque like the Gilbert type is prerequisite for 
the steady domain wall motion under magnetic fields. 
Energy dissipation via damping uses up the Zeeman en- 
ergy and sets the DW in motion. On the other hand, the 
DW under spin currents absorbs spin angular momentum 
via the bound magnon from conduction electrons and can 
move even without damping. 

Our theory is based on the perfect FM nanowires with- 
out spin damping. With the fully efficient absorption of 
spin angular momentum, the spin current sets the DW 
in motion with the adiabatic velocity v s , the spin waves 
are decoupled from the DW motion, and no Doppler ef- 
fect is expected in the spin wave energy spectrum. Note 
that the spin wave energy shift under spin currents was 
observed [22 ] in FM nanowires with uniform magnetiza- 
tion. According to the phenomenological LLG equation 
0) EH , the DW velocity is modified from the adia- 
batic value v s by the so-called a [23| and 3 fl3l damping 
torques. In this case more careful study [24| is required 
for elucidation of the Doppler effect in the spin wave en- 
ergy spectrum. 

In summary we studied the domain wall motion under 
spin current, based on the s-d model Hamiltonian. We 
found that the bound magnon plays an important role 
in the domain wall dynamics. Since the bound magnon 
is localized to the domain wall and has zero excitation 
energy, a spin current, without energy cost, transfers spin 
angular momentum to the domain wall via the bound 
magnon. Furthermore the hard axis anisotropy confines 
the bound magnon on the easy plane so that local spins 
can absorb spin angular momentum from the spin current 
and rotate about the hard axis without tilting away from 
easy plane. The bound magnon acquires its inertia due 
to the hard axis anisotropy and remains a free particle 
even under spin currents. With the full transfer of spin 
angular momentum from the spin current to the DW, 
spin waves are decoupled from the DW motion and no 
Doppler shift is expected for the spin waves. The bound 
magnon-mediated spin transfer mechanism leads to the 
free motion of bound magnon or domain wall under spin 
currents. 
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